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We discuss particle production in the high energy limit of QCD. Due to a large gluon density, the interaction 
reaches the black disk limit and the projectile is resolved into its partonic structure at the saturation scale. This 
leads to suppression of forward particle production and hereby to a faster absorption of cosmic ray air showers. 
This property is most suitable for the distinction of evolution scenarios for the saturation scale, e.g. fixed and 
running coupling BFKL, the latter of which is favored by air shower measurements. 



1. Introduction 

The biggest uncertainty in air shower simula- 
tions is certainly the hadronic interaction model, 
since QCD is poorly understood at these high en- 
ergies and accelerator data is not available. 

Higher twist corrections become increasingly 
important at high energies. Attempts to account 
for this are the implementation of an energy- 
dependent p t cutoff P for hard scattering or the 
resummation of enhanced pomeron diagrams in 
an efficient manner This way one tries to 
extend the applicability of hadronic interaction 
models up to GZK energies, E ~ 10 11 GeV. 

Our approach 3 is to consider the black 
disk limit (or black body limit - BBL) at high 
gluon densities within the Color Glass Conden- 
sate (CGC) approach 4, where the interaction 
probability is close to unity. This leads to a sup- 
pression of forward particle scattering, the most 
important phase space region for cosmic ray air 
showers. 

2. Hadron nucleus scattering at very high 
energies 

The elastic and total scattering cross sections 
for quark-nucleus scattering may be written as 
follows 0: 

a cl = J 'd 2 6[l-cxp(-Q 2 /47rA 2 )] 2 (1) 

a tot = 2 J d 2 b[l-cM-Q 2 s /^A 2 )] . (2) 

If Q s is large, the cross section approaches the 
geometrical limit. 
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Figure 1. Schematic view of the break-up of a 
hadron: the partonic structure is resolved at the 
saturation scale Q s 



Here, two different evolution scenarios are con- 
sidered for the saturation scale Qs- fixed coupling 

Ql(x,A)=Q%(A)^y , (3) 

and running coupling [E] 

Q 2 S = A 2 exp(log(g 2 /A 2 ) v /l + 2 cas2/ ) (4) 

ftp 

as ~ log(QVA 2 ) ' (5) 

with y — log(l/x). We assume Qo(A) to be pro- 
portional to the number of participants N part . 
The constant c assures a smooth transition to the 
fixed coupling scenario at low y. At high energies 
Q s s=s 5 GeV for running coupling and 20 GeV for 
fixed coupling for a central p-N collision 3 . 
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Figure 2. Suppression of forward scattering for 
central p-N collisions. 



The diagrammatic structure of an BBL event 
is schematically given in Fig. ^ The projectile 
interacts with the target as a whole and loses its 
coherence |7j. The partons are resolved at a scale 
given by the saturation momentum. An impor- 
tant feature is that soft physics is mostly sup- 
pressed, the typical transverse momentum being 
the resolution scale Q s . 

3. Monte Carlo implementation 

Since the PDFs are resolved at the satura- 
tion scale Q s , they provide the probability dis- 
tribution for a parton to appear in the final 
state. The mean transverse momentum is approx- 
imately given by Q s : 

Pi{x) = fi(Q 2 8 (x),x) (6) 
(Pt) « Qs(x) (7) 

Note that Q s is a function of x and b. For va- 
lence quarks we use the GRV94 parton distribu- 
tion functions. For a baryon, first the three va- 
lence quarks are generated, then the remaining 
energy is used for gluon bremsstrahlung. Here 
we use the ansatz from Ref. f5| 

xg(x, q 2 ) ~ — min (q 2 t , Q 2 s (x)) (1 - x) 2 , (8) 
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Figure 3. Mean multiplicity of central p-N 
events. 



which exhibits no divergence for low q t . The pro- 
duced gluons are ordered in rapidity and placed 
on strings between the valence quarks and the 
target nucleus, whose precise configuration is not 
important, since we are interested in forward scat- 
tering. In principle, a baryon-nucleus collision 
produces 3 strings, a meson- nucleus collision 2 
strings. However, when the invariant mass be- 
tween two of the three quarks is small, one can- 
not assume anymore that they fragment indepen- 
dently Therefore, we implement a cut-off in in- 
variant mass, 

m cut =m p = 0.77 GeV (9) 

below which two leading quarks are allowed to 
form a diquark. This is an important feature since 
it recovers the leading particle effect for low Q s . 
The hadronization of the strings is done within 
the LUND model 0. 

Of course, even at high energies, not all mini- 
mum bias events will be near the BBL. Therefore, 
we couple our model to a standard pQCD leading 
twist event generator, SIBYLL 2.1 1 . Which of 
the two models is to be applied for a given impact 
parameter, follows from the condition 
gNucicus^ _ 001 ) > lQe y (1Q ) 

The saturation momentum of the target at a lon- 
gitudinal momentum fraction of x = 0.001 for the 
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Figure 4. Inelasticity as a function of energy. 



projectile must be greater than a given value. The 
typical resolution scale for the valence quarks is 
correspondingly higher, since they have x w 0.2. 

4. Results of the BBL event generator 

Using the BBL as event generator, we can com- 
pare results in particle production to the pQCD 
model SIBYLL. At the same time, we compare 
the two different evolution scenarios, fixed and 
running coupling. Fig. [3 shows a Feynman-x dis- 
tribution of charged particles for central proton- 
nitrogen collisions. One notices the absence of 
particle production in the very forward region, 
which is the typical property of this approach. 
This figure shows running coupling only; fixed 
coupling would be even more extreme. The ab- 
sence of forward scattering is compensated by a 
relatively large multiplicity in the mid-rapidity 
region. In Fig. [31 we show the average multiplic- 
ities for fixed and running coupling BBL, as well 
as SIBYLL and QGSJET01 [H] results. Due to 
a missing ad-hoc pt cutoff, QGSJET01 exhibits 
also a large multiplicity. 

An interesting observable for air showers is the 
inelasticity, 
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where E max is the energy of the most energetic 



Figure 5. The fraction of BBL-events as a func- 
tion of energy. 



particle. A large value K w 1 means that most 
of the energy is used for particle production. A 
small value K « is significant for elastic or 
diffr active events, where most of the energy re- 
mains in the leading particle. For air showers, a 
high inelasticity means that the shower are ab- 
sorbed more rapidly in the atmosphere, giving 
rise to smaller A max values. Fig. 0] shows how 
the suppression of the leading particle results in 
a higher K, especially for fixed coupling. 

Finally, we show the fraction of BBL-events as 
a function of primary energy in Fig. It reaches 
90% at GZK energies. Note, that this does not 
mean that 90% of the cross section can be approx- 
imated as a black disk, since this also depends on 
the Q s with which the event was created. But 
it gives some idea about how this phenomenon 
becomes more important at high energies. 

5. Application of BBL to air showers 

We implement the BBL event generator into 
the Seneca air shower code an d compute the 
mean A max as a function of energy for proton 
primaries. 

In Figure we compare the predictions of 
SIBYLL 2.1 for proton and iron induced show- 
ers to the saturation model (BBL, for proton 
primaries only) with running and fixed coupling 
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Figure 6. X max as a function of primary energy. 



BFKL evolution of Q Sl respectively, and to pre- 
liminary Hires stereo data ^2j- In the saturation 
limit, showers do not penetrate as deeply into the 
atmosphere. This is due to the "break-up" of the 
projectile's coherence together with the sup- 
pression of forward parton scattering (for central 
collisions). The comparison to the data might 
suggest a preferably light composition, but the 
uncertainties at these energies are still consider- 
able. 

6. Conclusions 

We developed an hadronic interaction model 
considering the black disk limit in high density 
QCD. The suppression of forward scattering re- 
sults in faster absorption of particles within air 
showers, which leads to a smaller X max . This 
feature allows one to distinguish between possible 
evolution scenarios for the saturation momentum. 
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